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We have obtained high resolution (<, 1 ") deep x-ray images of the globular 
cluster 47Tucanae (NGC 104) with the Chandra X-ray Observatory to study 
the population of compact binaries in the high stellar density core. A 70 kilosec 
exposure of the cluster reveals a centrally concentrated population of faint 
(L^ ~10 30-33 erg s _1 ) x-ray sources, with at least 108 located within the central 
2' x 2.5' and ,>half with L x ^lO 30 5 erg s _1 . All 15 millisecond pulsars (MSPs) 
recently located precisely by radio observations are identified, though two are 
unresolved by Chandra. The x-ray spectral and temporal characteristics, as 
well as initial optical identifications with the Hubble Space Telescope, suggest 
that >,50% are millisecond pulsars, ~30% are accreting white dwarfs, ~15% are 
main sequence binaries in flare outbursts and only 2 to 3 are quiescent low mass 
x-ray binaries containing neutron stars, the conventional progenitors of MSPs. 
An approximate upper limit of ~470M Q for the mass of an accreting central 
black hole in the cluster is derived. These observations provide the first x-ray 
"color-magnitude" diagram for a globular cluster and census of its compact 
object and binary population. 

As the oldest stellar systems in the Milky Way Galaxy, globular clusters are laboratories 
for studies of stellar and dynamical evolution. The advantages of measurable ages and 
distances, both of which can be determined currently to ~10%, make globulars an especially 
attractive site for study of the evolution of stellar populations. The high stellar densities 
found in the cores of many globulars, with central values as high as 10 6 M Q pc~ 3 (1), 
implies that stellar interactions or near encounters can be relatively common - especially 
between stars and binary star systems. The role of binaries in cluster cores is dramatic: as 
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a source of localized gravitational binding energy, they act as a dynamical heat source in 
the cluster core, stabilizing it against further collapse into a central black hole (2). 

Compact binary stars, in which one member is either a white dwarf (WD) or a neutron 
star (NS), are markers of the extremes of stellar and binary evolution. Together with 
short-period binaries containing main sequence stars, compact binaries interact with stars 
in globular cluster cores and affect the stellar populations directly. Stellar evolution within 
compact binaries containing either a WD or NS leads to mass transfer, and accretion onto 
the compact object makes them (and their progeny, MSPs) visible as persistent x-ray 
sources. X-ray flare emission, and fainter continuous emission, can also be detected from 
chromospherically active main sequence (M-S) binaries, known as BY Draconis systems. 
Thus x-ray observations are a particularly efficient way to study the compact binary 
population and the WD vs. NS population in globulars. 

A complete sample of luminous x-ray sources (L x (2-10keV) >,10 36 erg s -1 ) has been 
discovered in globular clusters in the Galaxy: one bright source is detected (though nearly 
50% are transient) in each of 12 clusters, most with high central stellar density and 
metallicity (3). Type I x-ray bursts (4) have been detected from each, indicating all are NSs 
accreting from lower mass binary companions. This was also suggested by the ~3 "-accuracy 
x-ray positions measured for 6 of them with the Einstein x-ray Observatory which led to 
subsequent optical identifications for several and which yielded radial offsets consistent with 
their being ~2M Q objects (5). These 12 low mass x-ray binaries (LMXBs) in the system 
of ~150 globulars in the Galaxy are, as originally suspected (6), ~200x more common 
(per unit mass) in globulars than in the Galaxy as a whole. Thus 2-body stellar (tidal) 
interactions between a NS and M-S star, and the more probable exchange interactions of a 
NS with a pre-existing binary of two M-S stars, enhances LMXB production in globulars. 
The likely descendants of LMXBs, MSPs, in which the NS has been spun up to millisecond 
spin periods by accretion and becomes visible as a radio pulsar after the accretion phase 
ends (7), are also favored in globulars. The relative LMXB vs. MSP populations in a 
globular can thus constrain the formation epoch of each and whether alternative production 
channels are required for MSPs, such as accretion induced collapse (AIC) of massive WDs 
(to NSs) which could directly produce MSPs without a progenitor LMXB phase (7). In the 
relatively massive and high central density cluster 47Tuc (NGC 104), at least 20 MSPs are 
now detected which imply a total population of ~100-200 in this cluster alone (8). 

We obtained a 70 ksec exposure of 47Tuc on UT16.31 - 17.22 March, 2000, with the 
Chandra X-ray Observatory and Advanced CCD Imaging Spectrometer-Imager (ACIS-I) at 
its focus (9). Our primary goal was to study the population of low luminosity x-ray sources 
(L x <40 33 ^ 34 erg s -1 ) discovered in this and other globulars (10) and originally suggested 
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to be primarily accreting WDs and thus cataclysmic variables (CVs), together with a 
smaller population of the typically transient NS systems (LMXBs) in quiescence (quiescent 
LMXBs). Alternatively, the quiescent LMXB population may dominate the faint x-ray 
source population in globulars (11), and be required to maintain the large MSP population 
if they are derived continuously from LMXBs (7). Thus measurement and identification 
of a large sample of low luminosity x-ray sources in 47 Tuc, with its rich MSP population, 
can constrain the formation and evolution of MSPs as well as the relative populations of 
cluster WDs vs. NSs, with implications for the cluster initial mass function. We present 
results from a 2.0' x 2.5' field centered on the cluster core and chosen to include all of the 
15 MSPs with precise positions (12). This analysis for an inscribed radius of 1' = ~2.6r c , 
where r c is the cluster core radius (13), should include most compact binaries given their 
expected mass segregation (2). Results for the full 16' x 16' ACIS-I field will be presented 
subsequently. 

Overview of Source Distributions 

Previous x-ray studies of 47Tuc with the ROSAT telescope and HRI detector (14) 
revealed a population of 9 sources with x-ray luminosity ^3-100 x 10 erg s within 
~1 ' of the cluster center as well as underlying unresolved emission with total luminosity 
L x ~4 x 10 32 erg s -1 . The factor of >10 improvements in angular and spectral resolution 
(to ~0.5" and ~120eV, respectively, with ACIS) achieved by Chandra over ROSAT enable 
detection of sources ~10-100x fainter. The spectral resolution enables comparison with 
x-ray spectra of known classes of source, while the exquisite positional resolution permits 
much more detailed searches for optical counterparts with high resolution HST images. Our 
deep exposure reached a sensitivity limit of L x ~6 x 10 29 erg s _1 for a threshold source with 
3 counts (vs. ~0.5 counts total background per source) detected at 90% confidence by the 
Chandra analysis tool wavdetect (15) and assuming a soft thermal bremsstrahlung spectrum 
with temperature kT = 1 keV (see below) for comparison with the ROSAT results (14). 
The "true" color image of the cluster core as imaged by Chandra (Fig. 1) shows the range 
of source brightness and spectral color by combining counts detected in soft (0.5 - 1.2 keV; 
red), intermediate (1.2 - 2 keV; green) and hard (2-6 keV; blue) bands. Wavdetect finds 
103 sources with total counts (in a 0.5 - 4.5 keV primary detection band) from 3 (detection 
threshold) to 5332 over the total 70 ksec exposure. Five additional sources are found at 
the positions of MSPs Q, L, T, M and C: Q, L, and T were missed by wavdetect due to 
crowding, while M and C are only 3 and 1 count "detections" , though precisely at the MSP 
positions, due to low exposure (~30%) on the detector gaps. Using the recently derived 
distance (5.0kpc) and absorption for 47Tuc (16), we derive the luminosity distribution for 
the sources (Fig. 2). The complete list of 108 source positions and counts detected in 
several bands, as well as source identifications currently derived (some discussed below), are 
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given in the Supplementary Table (available at Science On Line). 

The x-ray luminosity function increases at low values, with more than half the sources 
with L x <3 x 10 30 erg s _1 . The 23 sources in the lowest L x bin (2.5 - 7.5 cts; L x ~10 30 
erg s" 1 ) of the linear distribution form an incomplete set. From visual examination of 
the image as well as a ~500 count excess of total counts over detected source counts in 
the central core region (r <lr c ), we estimate an additional >,100 sources with L x <10 30 
erg s _1 are missed with wavdetect due to crowding. Both a statistics-limited overlapping 
source detection algorithm (17) and a deeper ACIS-S exposure (with better low energy 
sensitivity) are needed to constrain the number and types of these faintest core sources. 

The sources appear to be of (at least) 4 types as marked with different symbols in Fig. 
1 and discussed below. The distribution of source types can also be seen in the x-ray "color 
magnitude" diagram (Fig. 3), which is possible to derive for the first time for a significant 
population of x-ray sources in a single globular cluster. 

Source Identifications 

Our analysis region of the cluster core was chosen to include the 15 MSPs with known 
positions (12) because 10 MSPs (including one in the globular cluster M28) have been 
detected in x-rays (18) and detection of even some in 47Tuc would enable a study of a 
population of cluster MSPs and their relation to both quiescent LMXBs and CVs. All 
14 of the resolvable MSPs are detected (MSPs G and I, possibly a bound triple, have 
only 0.12" separation (12)), with exposure-corrected counts ranging from 3-27 and thus 
L x ~10 29 - 8-30 - 8 erg s" 1 . Source identifications are possible, since even the faintest (3ct; >J2a 
in wavdetect) Chandra sources have positions with ~0.2" uncertainties (la) due to the 
~0.8" Chandra image diameter (near the telescope axis) and 0.5" ACIS pixel size. The 
precise (~0.001") radio timing positions of the MSPs permit an astrometric solution for the 
Chandra source positions. Using 7 of the brighter or best resolved Chandra/MSP source 
candidates, the rms deviation between the Chandra and precise MSP positions is only 0.11 ". 
Given a similarly precise Chandra-HST astrometric solution from our optically identified 
sources (see below), this will permit optical searches for the MSP counterparts. Because 
the 14 resolved Chandra MSP candidates are all located near the center of the circles which 
mark the precise MSP positions (Fig. 1), the identifications are very likely. Although the 
probability that one of the estimated ~100 unresolved sources in the 23" radius core is 
within 0.5 "of any MSP position is 0.05, this is applicable to only MSPs L and T (the only 
MSPs in the core not detected automatically with wavdetect, and thus with positional 
uncertainties ~1 pixel = 0.5" instead of wavdetect centroid values of typically <,0.2"). 
The ACIS time resolution (nominal 3.2sec) does not permit analysis for pulsations, and the 
limited counts do not permit spectral analysis other than the determination of hardness 
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ratios. The relatively narrow x-ray luminosity range implies a steep x-ray luminosity 
function, with typical L x similar to the field MSPs detected outside of clusters (18), though 
the predominantly "red" x-ray colors (see below) are in contrast to some of the field MSPs 
with hard spectral components and the one in M28 with L x ~10 32_33 erg s _1 (18). Full 
details of the 47 Tuc MSP x-ray source characteristics will be presented separately (19). 

Three other classes of sources are suggested by their x-ray spectra and temporal 
characteristics (Fig. 1). X5 and X7 are probable quiescent LMXBs, 13 sources are candidate 
CVs, and 6 sources appear to be main sequence (M-S) binaries detected in flare outbursts 
due to enhanced chromospheric activity. Each class is considered in turn to then enable 
constraints on the much larger unidentified population. 

X5 and X7 are detected in the 0.5-4.5keV medium-energy band with 4435 and 5332 
counts (hereafter medcts) and have spectra dominated by soft blackbody components with 
kT= 0.31 and 0.29 keV, respectively. Their apparently constant black body luminosities 
(log L x = 32.8; comparable to that detected with ROSAT (14)) and implied emission radii of 
only ~0.8 km suggest they are quiescent LMXBs in which the x-ray emission is dominated 
by incandescence of the hot NS, slowly cooling from a previous accretion outburst phase. 
Using a hydrogen-atmosphere model, which modifies the opacity and temperature, as 
demonstrated with Chandra ACIS-S spectra for the field quiescent LMXB Cen X-4 (20), 
more realistic NS radii of ~10 km are derived (21). Our detailed spectral analysis of X5 
and X7 also reveal a power law component (photon index 2.6 - 3) similar to that found for 
Cen X-4 and indicative of either residual accretion or a "propeller" spin-down luminosity 
source, as well as a possible line emission component (modeled (21) as a collisionally excited 
Raymond-Smith plasma with kT = IkeV), indicative of a thin corona or perhaps wind from 
the system. The integrated fluxes for these two components are ~0.3 to 0.5 and ~0.2 of 
the black body component, respectively. Although a statistical test (Kolmogorov-Smirnov) 
on photon arrival times of X7 shows no significant variations in flux, X5 shows dramatic 
dips and possible eclipses, which may also be detected in the optical counterpart we have 
identified with HST (22). Both X5 and X7 are most probably quiescent LMXBs given 
their similar luminosities and spectra to those tabulated (23) for the 6 known NS-quiescent 
LMXB systems in the field. However, the hot gas and possible wind component may 
indicate that X5 and X7 are luminous newly-born MSPs ablating their companions and still 
optically thick to their radio emission, because they are not among the currently located 
MSP sample (12). The lack of positive flux increases (i.e. flares; vs. the negative dips seen 
in X5) and otherwise constant flux would support, but not require, the MSP interpretation. 

The 13 candidate CVs (cf. Fig. 1) identified thus far in our HST identification program 
(22) have blue stellar counterparts indicative of accretion disks (though not all are yet 
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determined to be variable). Apart from the quiescent LMXBs X5 and X7, all remaining 
ROSAT sources [X6, X9, X10, Xll, X13, and X19; (14)] within our analysis region (X4 is 
just outside) are CV candidates. Spectra have been fit (21) for the brightest 8 sources (with 
medcts = 135 - 2402) and most are best fit with thermal bremsstrahlung with temperatures 
kT = 5-30 keV, typical of CVs (24). Soft black body components, as found for quiescent 
LMXBs, are generally not present. Three sources have unusually hard spectra: X6 is best fit 
with a thermal bremsstrahlung spectrum with kT >,75 keV, and both W8 and W15 require 
similarly hard spectra but with internal self-absorption column densities of NH = 1.4±0.4 
and 0.7±0.2 x 10 22 cm' 2 , respectively. Power law fits are equally acceptable for all three 
but with photon index a ~1.1±0.1 (consistent with the slope of a thermal bremsstrahlung 
spectrum). Since such a flat power law index at low energies (<6 keV) is unprecedented 
for MSPs (for which the limited spectral data (18) indicate a ~2 - 3) or active galactic 
nuclei (AGN) with typical (25) a ~1.7, a thermal bremsstrahlung model is more likely. The 
excess NH in sources W8 and W15 is evident from their "blue" colors in Fig. 1 (absorption 
of soft counts) and might suggest they are background AGN, probably Seyfert 2 galaxies, 
even though the number expected (from hard sources in Chandra deep surveys (26)) in this 
2' x 2.5' field and exposure is <,0.1. However, not only are their power law spectral fits too 
hard, but our identification (22) of both W8 and W15 with blue stars imply x-ray/optical 
flux ratios similar to values expected for CVs (e.g. \og{¥ x /Y v ) = +0.1 for W15). X-ray 
self-absorption is likely in high inclination CVs (27) and may also arise from the "accretion 
curtain" (27), or relatively cool inner disk near the magnetospheric-disk interface suspected 
to exist in magnetic CVs in which accretion onto the WD is channeled by the magnetic 
field. Magnetic CVs also typically have relatively high temperature thermal bremsstrahlung 
components, though the kT values for X6, in particular, is unusually large. 

The brightest CV candidate, X9 (medcts = 2402), is well fit with a thermal 
bremsstrahlung spectrum with kT = 32±6 keV and is coincident (0.1") with the blue CV 
candidate VI (28), confirming the identification by the ROSAT detection (14) of X9. The 
70 ksec Chandra observation of X9 shows evidence for pulsations (false alarm probability 
~2 x 10~ 3 ) with period P = 218.24s and modulation depth 30% (Fig. 4), suggesting this 
is direct evidence for an accreting magnetic WD and thus magnetic CV. Although the 
total flux from X9 is relatively constant (Fig. 5), and a factor of ~2 fainter than detected 
with ROSAT, the pulsation amplitude is variable (strongest in the second quarter of the 
observation) and could instead indicate transient pulsations as found in some CVs (24). The 
faint CV candidate AK09 is also a possible magnetic CV. It was not detected with ROSAT 
and is too faint even with Chandra (medcts = 36) for spectral fitting but is relatively soft 
(red; cf. Fig. 1). The Chandra position is coincident (0.1") with the blue variable AK09 
(29) with l.ld binary period, and originally suggested (30) as the counterpart of the single 
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x-ray source detected in the original Einstein x-ray image of 47 Tuc (5), for which the 
positional uncertainty (3"; 90%) included both X9 and AK09. HST spectra (31) now show 
AK09 to be almost identical to the magnetic CV system GK Per, with similarly long (2d) 
binary period and occasional outbursts. AK09 also shows UV outbursts (32) so that the 
bright Einstein source, which had both a thermal bremsstrahlung spectrum (kT = 2-6±l 
keV) as well as possible 120.2s pulsations (30), might have been AK09 in outburst. Both 
Chandra (21) and HST (22) studies further constrain the magnetic CV hypothesis. 

The optically discovered CV, V2 (33), which displays dwarf nova-like outbursts is 
identified (0.1") with source X19 and fit with a thermal bremsstrahlung spectrum with 
kT = 12 ±6 keV, and the blue variable V3 (32) is the probable counterpart (0.1") for 
X10, in contrast to the claim from ROSAT data (14). The Chandra spectrum of X10 is 
unique among all the (moderately) bright sources for which spectral fits were possible. X10 
can be fit with a pure power law (a = 2.7±0.2). Its x-ray color and extreme variability, 
with dips or flares and eclipses (Fig. 5), for which a power spectrum reveals a significant 
3.83h period, are similar to the quiescent LMXB candidate X5 although a black body + 
power law component (as for X5) is only acceptable if a Raymond-Smith component is 
added in absorption. However, its blue variable counterpart V3, with Balmer emission 
lines (34), suggests accretion is ongoing. We suggest X10, and thus perhaps X5, may be 
similar to the peculiar magnetic CV AE Aqr, in which spin energy of the WD is partly 
converted to non-thermal radiation by a magnetic propeller, possibly anchored in the disk 

(35) . Alternatively, the power law spectrum suggests X10 might be, as mentioned for X5, 
an enshrouded MSP. 

The fourth group of probable source identifications are main sequence binaries. Of the 
6 sources marked (Fig. 1), the brightest, identified with the 12. 7h binary E8 (29), is the 
most convincing association with what may be chromospherically active M-S binaries. No 
blue excess is found in our HST analysis (22), and the star appears just above the main 
sequence. This is as expected for a BY Dra star, the main sequence analog of RS CVn 
systems in which x-ray emission is produced primarily in giant flare outbursts. Although 
the L x = 10 31 ' 4 erg s _1 is a factor of ~10 larger than the maximum found with ROSAT 

(36) for field BY Dra systems in quiescence, this is probably due to what appears to be 
a smooth decline from a large outburst (Fig. 5). The other 5 candidate BY Dra systems 
(Fig. 1) also show x-ray flare-like variations. Using our Chandra-HST astrometry, we have 
found at least 10 other optical variables without blue excesses that are possible BY Dra 
systems, although without x-ray flare variations they are less likely identifications. Few, if 
any, RS CVn systems are seen, with only 2-3 possible matches (<,0.5") of Chandra sources 
with sub-giants. 
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Implications for Population of Compact Objects, Binaries and Cluster Core 

While source types for some individual cases are uncertain (e.g. X10), it appears 
that the 108 sources near the core of 47Tuc are a mixture of 4 source types. Their 
quasi-continuous distribution in x-ray colors (Fig. 3) is due to the relative insensitivity 
of broad-band colors to detailed spectral differences, particularly for multi-component 
spectra. Nevertheless, some broad differences emerge when the distribution of x-ray color 
is examined (Fig. 6). Most striking is that most of the MSPs (for which identifications 
are unambiguous) are relatively soft: 8 of the 12 with medcts >,5 (and thus rough colors 
available, as in Fig. 3) are softer than any of the much brighter CV candidates. Similarly, 
the two most probable quiescent LMXBs are softer than all of the CVs (except AK09) 
and are intermediate between the CVs and MSPs. The BY Dra candidates (BYs) are 
intermediate in color but closer to the CVs. Although the statistics are limited, the values 
for the mean, Xcolor , and standard deviation, a xco i or , for each population are: (Xcolor , 
(Jxcoior ) = (1.8±0.3, 0.98), (-0.1±0.3, 1.09), (0.8±0.2, 0.56) and (1.5±0.1, 0.17) for MSPs, 
CVs, BYs, and quiescent LMXBs respectively. The corresponding distribution values for 
the unidentified (UnID) sources are (1.1±0.1, 1.03), or between the BYs and MSPs. The 
BY contribution to the UnlDs can be measured by a deep followup observation. Most of the 
same UnID sources will not be detected if they are flare outbursts from BY Dra systems, 
while MSPs will be constant in their x-ray emission. Given the approximate relative 
numbers of the identified sources and the Xcolor distributions (Fig. 6), we decompose 
the 61 UnID (>5ct)source distribution into 2 quiescent LMXBs, 18 CVs, 11 BYs, and 30 
MSPs. Combined with the identified sources, this yields approximate relative (fractional) 
contributions of 0.04, 0.3, 0.16 and 0.5 for the populations of quiescent LMXBs, CVs, BYs, 
and MSPs with L x >,10 30 erg s" 1 , respectively, where the uncertainties are probably 30%. 

We conclude the source population in 47Tuc is dominated by MSPs, with a total 
number of >, 45-60 in the present Chandra sample (with L x >1 x 10 30 erg s _1 ). These 
would naturally account for the background or unresolved red sources, suggesting a still 
larger MSP population at lower luminosities, if these faint sources are detected as constant 
flux sources in followup observations. Further inferences for the MSP population, including 
their spindown energy loss, E, vs. L x correlations and radial distributions in the cluster, 
are discussed separately (19). However, the CV population may also extend into this faint, 
soft source background population if there are large numbers of strongly magnetic (AM Her 
type) CVs for which soft black body emission dominates. These would be optically very 
faint (like the MSPs), but would display marked variability (relatively long high- low states), 
and would thus also be distinguished from both BY Dra and MSP systems of comparable 
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Our estimate that the present sample contains ~30 CVs is ~l/3 the number predicted 
for tidal capture CV production models (37). This may indicate tidal capture is relatively 
inefficient, but in fact the prediction of a large excess of low luminosity sources (L x <,10 30 
erg s _1 ) may also be consistent with our "background" sources. The size of the CV sample 
is unexpected because so few dwarf novae have been found (33). Combined with the 
evidence presented here that the brightest CV candidates (e.g. X9) are magnetic CVs, it 
may be that CVs in globulars are indeed dominated by magnetic systems, as suggested 
from HST spectra (38), in which dwarf nova outbursts are suppressed, at least for short 
binary period systems. 

If only 2-3 quiescent LMXBs (X5, X7, and possibly X10) are present in the core with 
the large sample of ~100 MSPs inferred here or from radio counts (8), the spin-up timescales 
for MSPs are reduced to T sp inu P ~N quiescentLM xBs/N M sPs x Tmsps, where Tmsps ~3 x 10 9 
years is the typical MSP spindown age. Thus T sp j„ up ~1 x 10 s years, which is ~10x faster 
than the >,10 9 years needed to transfer the ~0.03P4~ 4 ^ 3 M Q required (7) to spin up the 
NSs in 47Tuc to their typical ~4msec spin period (P4 is the final period in units of 4msec) 
at the time-averaged mass accretion rate fn ~3 x 1O _11 M yr _1 appropriate to typical 
LMXBs in clusters with L x ~K) 36 - 3 erg s _1 (3) which may be in quiescence ~50% of the 
time. Thus the observed quiescent LMXBs are insufficient to support the MSP population 
and alternative paths to MSP production may be required such as direct collapse of massive 
WDs in high m CVs (7) (for which X9 may be relevant) or MSP formation from common 
envelope evolution of a NS with a prior generation of ~1-3M Q binary companions leading 
to a short-lived NS-WD LMXB phase (39). 

Finally, the Chandra image allows the most sensitive mass limits for a black hole in 
the core of a globular cluster. The precise MSP positions and associated pulsar dispersion 
measures (12) imply an electron density (and thus ionized gas at T ~10 4 K) of n ~0.1 
cm -3 in the cluster core. Bondi-Hoyle spherical accretion with radiative efficiency e onto a 
compact object of mass M, radius R and velocity V in a gas with number density n and 
sound speed c s gives an accretion luminosity L acc = e G M 3 m p n (V 2 + c 2 )~ 3 / 2 R _1 , where 
m p is the proton mass. For an isolated NS as recently considered for low luminosity x-ray 
sources in globulars (40), this gives L acc = 10 30 erg s~ 1 for 47Tuc and appropriate NS values 
of M ns = IAMq , R ns = 10 km, e = 0.1 and V ns ~c s ~10km/s. If the accretion energy 
is thermalized over the NS surface, the expected black body spectrum with kT ~15 eV is 
too soft to account for the detected Chandra sources. Although the thermalization process 
is uncertain, higher temperature radiation from a magnetic pole is also inhibited at such 
low accretion rates in the spherical case if the NS is rapidly rotating. However, for a black 
hole, with larger mass M^, emission radius ~3 Schwarzshild radii and lower equipartition 
velocity V^ = V„ s (M ns /Mf ) / l ) 1 / 2 (so that c s dominates), the accretion luminosity will 
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likely arise from a much lower efficiency flow with e ~10 -4 and radiate a much higher 
temperature spectrum (kT ~100 keV), such as calculated in advection (or convection) 
dominated accretion models (41) for accretion on the black hole x-ray source SgrA* in the 
galactic center. We estimate a soft x-ray luminosity L x (0.5-2.5keV) ~4.5 x 10 25 e_4M| /l Tioo 
erg s _1 , where the uncertain radiation efficiency e and spectral temperature T are scaled 
to 10 -4 and 100 keV, respectively. The Chandra upper limit of L x ~1 x 10 31 erg s _1 for 
the brightest source in the cluster center error circle (42) (Fig. 1) implies an upper limit of 

— 1/2 

Mbh ~47OM e__ 4 Tioo- This is consistent with the understanding that binary heating (2) 
prevents core collapse into a moderately massive central black hole and that stellar mass 
(~3-10 M Q ) black holes produced from the initial stars more massive than the NS cutoff 
have not coalesced in the core but rather have been ejected by the hardest binaries (43). A 
more precise cluster center and accretion model can improve these limits for M^, which 
are already below the 1700M Q value derived from the central surface brightness profile 
(44). Similar limits for those post core collapse clusters which, like 47Tuc, do not contain a 
bright LMXB and for which central gas densities could be estimated from MSPs, would be 
especially interesting. 
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Fig. 1. X-ray color image of the central 2' x 2.5' of 47Tuc. Colors are derived from counts 
detected in soft (0.5 - 1.2 keV; red), intermediate (1.2 - 2 keV; green) and hard (2-6 keV; blue) 
bands with the 13 chip of the ACIS-I detector on Chandra. The pixel size (visible in enlargement 
of central core) is 0.4914", and the enlarged central region is 35 "square. Source identifications 
are: MSPs (circles), marked with source letter names (12); quiescent LMXBs X5 and X7; CV 
candidates (squares), marked with ROSAT Xnumbers (14) or other Chandra Wnumbers; and 
possible flaring BY Dra systems, or M-S binaries (triangles), marked with E8 (29) or Chandra 
Wavdetect source numbers (Wnumbers). Precise positions, fluxes, hardness ratios, bright source 
spectra and variability are tabulated for all 108 sources shown here (73 unidentified are not 
labeled, for clarity) in our full-field analysis (21), with only key values given here in the text. 
The mean cluster center (41) and 3a error (3" radius dashed circle) are marked; other symbols 
(squares, etc.) are centered on sources but with sizes much greater than positional uncertainties. 
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Fig. 2. X-ray luminosity and counts distributions of sources detected with >2.5cts in central 
2' x 2.5' of 47Tuc. Counts detected in medium band (0.5 - 4.5 keV) are corrected for instrument 
and telescope response, 5.0 kpc cluster distance (16) and known interstellar absorption of 
Ay = 0.12 (or equivalent column density NH = 2.4 x 10 20 cm~ 2 ) and an assumed thermal 
bremsstrahlung spectrum with kT=l keV to give source luminosity in the 0.5-2.5 keV band for 
direct comparison with ROSAT results (14). The extended energy band of Chandra (0.5 - 8 keV) 
is utilized for hardness ratio and spectral analysis, which also enable conversion of L x between 
the Chandra and ROSAT bands. (A) log source counts vs. log L x , with log L x value given at 
the center of 0.3dex bins, (B) linear source counts vs. actual counts detected , with value given 
at center of 5ct wide bins. Corresponding log L x values range from 29.9 to 31.4. The 96 sources 
plotted in (B) account for only 13% of the total detected x-ray flux from the cluster central region, 
whereas the brightest two sources (probable quiescent LMXBs X5 and X7) contribute 57% of the 
total counts. 
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Fig. 3. X-ray brightness (medcts= 0.5-4.5keV) vs. color (ratio of soft (0.5-1. 5keV) /hard (1.5-6keV) 
counts) distribution for sources with medcts >5 cts. Source types are labeled for the sources 
identified with 12 MSPs (red •), 2 quiescent LMXBs (orange X), 13 CVs (blue A) and 6 
M-S binaries (BY Dra systems) in flare outbursts (yellow O). The 61 remaining (green *) are 
unidentified. Source hardness increases to left. 
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Fig. 4. Folded light curve for the brightest CV candidate, X9, showing its probable detection 
as a 218.24s pulsed source and thus identification as an accreting magnetic WD, or magnetic CV 
system (a NS x-ray pulsar with such a long period would be expected only in a wind-fed high 
mass x-ray binary found only near regions of massive star formation in the Galaxy, not in globular 
clusters) . 
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Fig. 5. Light curves, over the full 74ks (duration) observation, for 3 of the 4 source types (all 
MSP fluxes are constant, within statistics). (A) quiescent LMXB source X7; (B) CV candidates 
X9 (cf. Fig. 4), showing flickering, and (C) X10, showing dips and eclipses with 3.8h period 
(similar to X5); and (D) M-S binary candidate E8, showing smooth decline from probable large 
flare. 
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Fig. 6. Distributions of x-ray color (cf. Fig. 3) for 94 sources with medcts >,5 counts and 
proposed source identifications (see Figs. 1 and 3; source types are color coded same as in Fig. 3). 
M-S binaries are BY Dra systems in flaring outbursts. Means and standard deviations for each 
source distribution are given in the text. 



- 1 - 



Table 1. 47 Tuc Core X-Ray Sources (Grindlay et al 2001, Science On Line, 17 May 2001) 
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Table 1 — Continued 
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14.8 


5 


2.9 


4.4e30 




W62 


0:23:59.434±0.02 


-72:04:47.79±0.2 


9.5 


7.9 


2 


4 


2.0e30 




W63 


0:23:59.392±0.03 


-72:03:58.81±0.08 


10 


5.9 


5.9 


1 


3.0e30 


MSP-J 


W64 


0:23:57.699±0.007 


-72:05:01. 93±0.03 


164.2 


88.4 


87.4 


1 


3.4e31 




W65 


0:23:53.548±0.02 


-72:05:11. 20±0.1 


9.7 


6.8 


3 


2.3 


4.0e30 




W66 


0:24:16.023±0.03 


-72:04:51. 56±0.2 


10.3 


10.3 


0.5 


20 


2.2e30 




W67 


0:24:13.818±0.03 


-72:04:43.83±0.1 


13.9 


8.9 


5.6 


1.6 


2.9e30 


MSP-D 


W68 


0:24:13.477±0.02 


-72:05:23.42±0.2 


6.4 


5.4 


1 


5.3 


1.3e30 


M-S 


W69 


0:24:12.643±0.02 


-72:04:22.38±0.2 


3.9 


2.8 


2 


1.4 


8.2e29 




W70 


0:24:11.915±0.02 


-72:04:43.84±0.1 


11 


8 


3 


2.6 


2.3e30 




W71 


0:24:10.642±0.01 


-72:05:16.37±0.2 


8.6 


6.9 


1 


7.1 


1.8e30 




W72 


0:24:10.484±0.02 


-72:04:59.92±0.1 


16.6 


15.1 


3 


5 


3.5e30 




W73 


0:24:07.224±0.02 


-72:04:57.39±0.1 


13.5 


8.1 


2 


4 


2.8e30 


- 


W74 


0:24:06.686±0.02 


-72:04:06.74±0.1 


10.4 


10.6 


0.5 


21 


2.2e30 


MSP-H 


W75 


0:24:06.429±0.02 


-72:04:52.79±0.1 


24.6 


15 


8 


2 


5.2e30 




W76 


0:24:6.436±0.04 


-72:04:30.11±0.1 


5.1 


2 


4 


0.5 


l.le30 




W77 


0:24:03.92±0.02 


-72:04:42.71±0.08 


27.3 


19.8 


6 


3.3 


5.7e30 


MSP-F,S 


W78 


0:24:03.516±0.02 


-72:05:22.54±0.1 


4.2 


4.3 


0.5 


8.6 


8.8e29 




W79 


0:24:03.565±0.02 


-72:04:41. 50±0.1 


19.1 


14 


8.6 


1.6 


4.0e30 




W80 


0:24:02.658±0.02 


-72:04:41.47±0.1 


11.4 


8.5 


2 


4.2 


2.4e30 
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Table 1 — Continued 



Src a 


RA 6 


Dec b 


Medcts c 


Sfcts 


Hdcts 


Sf/Hd d 


L x e ID f 


W81 


0:24:01. 775±0. 03 


-72:05:12.30±0.1 


7.9 


4 


5.2 


0.8 


1.7e30 


W82 


0:24:01.434±0.02 


-72:04:41. 72±0.1 


8.6 


8.9 


0.5 


18 


1.8e30 


W83 


0:24:00.977±0.03 


-72:04:13.70±0.1 


9.5 


7.6 


3 


2.6 


2.0e30 


W84 


0:24:00.424±0.04 


-72:05:46.55±0.2 


6 


2 


3 


0.7 


1.3e30 


W85 


0:23:59.404±0.04 


-72:04:38.46±0.1 


9.5 


8.5 


1 


8.3 


2.0e30 


W86 


0:23:58.988±0.02 


-72:05:09.70±0.1 


3.8 


3 


1 


3 


8e29 


W87 


0:23:57.884±0.04 


-72:04:52.76±0.08 


7.4 


4.5 


3 


1.5 


1.6e30 


W88 


0:23:54.175±0.03 


-72:05:14.50±0.1 


4.6 


3 


2 


1.5 


1.7e30 


W89 


0:23:52.544±0.03 


-72:05:32.99±0.1 


3 


3 


0.5 


6 


1.6e30 


W90 


0:23:50.806±0.03 


-72:05:41. 51±0.1 


8.3 


5.4 


3 


1.8 


2.1e30 


W91 


0:24:17.55±0.05 


-72:03:58.28±0.1 


3.7 


0.5 


4 


0.13 


7.8e29 


W92 


0:24:15.678±0.02 


-72:04:41. 31±0.2 


9.2 


7 


3 


2.3 


1.9e30 


W93 


0:24:12.017±0.02 


-72:05:08.06±0.2 


3.8 


2 


2 


1 


8.0e29 


W94 


0:24:11. 78±0.03 


-72:05:14.61±0.09 


5.5 


4 


2 


2 


1.2e30 


W95 


0:24:10.861± 


-72:04:15.84±0.1 


2.9 


3 


0.5 


6 


6.1e29 


W96 


0:24:07.281±0.02 


-72:05:19.58±0.2 


4.4 


4 


1 


4 


9.2e29 


W97 


0:24:05.379±0.02 


-72:05:11. 12±0.1 


5.5 


5 


1 


5 


1.2e30 


W98 


0:24:05.056±0.03 


-72:05:02.77±0.1 


10.5 


11 


0.5 


22 


2.2e30 


W99 


0:24:02.37±0.04 


-72:05:50.66±0.09 


4.7 


4 


2 


2 


9.9e29 


W100 


0:24:01. 368±0.02 


-72:04:27.26±0.2 


5.4 


5 


0.5 


10 


l.le30 
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Table 1 — Continued 



Src a 


RA b 


Dec 6 


Medcts c 


Sfcts 


Hdcts 


Sf/Hd d 


L x e 




W101 


0:23:59. 974±0. 02 


-72:05:31. 93±0. 2 


5.4 


5 


0.5 


10 


l.le30 




W102 


0:23:57.687±0.05 


-72:05:42. 98±0.1 


2.8 


3 


0.5 


6 


5.9e29 




wi rn 


u. z,o .o i . yyuxu . 


79-n4-9Q 31 +f) 1 

1 Zi.Ut.ZiC/.OIXU. 1 


1 Q 


o 


yj.o 


A 


A flp9Q 




W104 


0:24:16.49 


-72:04:25.15 


8 


7 


1 


7 


1.7e30 


MSP-Q 


W105 


0:24:08.54 


-72:04:38.91 


6 


5 


1 


5 


1.3e30 


MSP-T 


W106 


0:24:03.77 


-72:04:56.90 


12 


11 


1 


11 


2.5e30 


MSP-L 


W107 


0:23:54.49 


-72:05:30.72 


3 


3 


0.5 


6 


2.3e30 


MSP-M 


W108 


0:23:50.35 


-72:04:31.49 


1 


1 


0.5 


2 


1.3e30 


MSP-C 



Note. — a Chandra source number. b RA (right ascension) and Dec (declination) coordinates are given 
as epoch J2000. Errors are formal la centroiding errors from Chandra x-ray detection routine wavdetect. 
The final five MSPs were not formally detected by wavdetect, so the radio positions of Freire et al, 
2001 (12) are given without errors. c Medcts=counts in 0.5-4.5 keV band, Sfcts=counts in 0.5-1.5 keV 
band, Hdcts=counts in 1.5-6 keV band. Values of 0.5 counts are estimated upper limits. d Softness ratio, 
defined as Sfcts/Hdcts [in Figs. 3 and 6, Xcolor = 2.51og(Sfcts/Hdcts)]. e Intrinsic luminosities in ROSAT 
(0.5-2.5 keV) band, with multiplicative exposure corrections (=1.4, 2.0, 1.7, 2.6, 1.2, 2.5 and 1.9 for 
sources W63(MSPJ), W65, W88, W89, W90, MSPC and MSPM, respectively), based on 1 keV thermal 
bremsstrahlung spectrum, N H = 2.4 x 10 20 cm~2, d=5.0 kpc. X number from ROSAT (14) observations 
(e.g. X5), certain optical or radio counterparts (e.g. VI, AK09, MSP-D), and class identifications: 
qLMXB=quiescent LMXB, M-S=main sequence binary, CV=cataclysmic variable. Detailed discussion of 
of individual x-ray source spectra and variability is deferred to Heinke et al. 2001 (21) and probable optical 
counterparts are discussed by Edmonds et al. 2001 (22). 



